Effet Réponse

Effet de la réponse dans toutes les cellules

1. Liste des gènes différentiellement exprimés

RCHOP

singlet <- NormalizeData(singlet, assay = "RNA")
singlet <- ScaleData(singlet, features = rownames(singlet))
response <- FindMarkers(singlet, ident.1 = item1, ident.2 = item2 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod")

response.RCHOP <- response
save(response.RCHOP, file = "../Reponse/response.RCHOP.raw.Rdata")

response <- SC.Result(response)

DT.table(response)
# ES <- enrichIt(obj = singlet, gene.sets = GS, groups = 1000, cores = 8, min.size = 5)
# save(ES, file = "../Reponse/RCHOP_ES.RDs")

Excipient

singlet.exp <- NormalizeData(singlet.exp, assay = "RNA")
singlet.exp <- ScaleData(singlet.exp, features = rownames(singlet.exp))
response.exp <- FindMarkers(singlet.exp, ident.1 = item1, ident.2 = item2 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod")


save(response.exp, file = "../Reponse/response.exp.raw.Rdata")

response.exp <- SC.Result(response.exp)

DT.table(response.exp)
# ES <- enrichIt(obj = singlet.exp, gene.sets = GS, groups = 1000, cores = 10, min.size = 5)
# save(ES, file = "../Reponse/Excipient_ES.RDs")

Pré-greffe

singlet.pg <- NormalizeData(singlet.pg, assay = "RNA")
singlet.pg <- ScaleData(singlet.pg, features = rownames(singlet.pg))
response.pg <- FindMarkers(singlet.pg, ident.1 = item1, ident.2 = item2 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod")

save(response.pg, file = "../Reponse/response.pg.raw.Rdata")

response.pg <- SC.Result(response.pg)




DT.table(response.pg)
ES <- enrichIt(obj = singlet.pg, gene.sets = GS, groups = 1000, cores = 10, min.size = 5)
## [1] "Using sets of 1000 cells. Running 21 times."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
## Setting parallel calculations through a SnowParam back-end
## with workers=10 and tasks=100.
## Estimating ssGSEA scores for 50 gene sets.
## [1] "Calculating ranks..."
## [1] "Calculating absolute values from ranks..."
save(ES, file = "../Reponse/Pre-greffe_ES.RDs")

2. Enrichissement

RCHOP

load(file = "../Reponse/RCHOP_ES.RDs")

singlet <- AddMetaData(singlet, ES)
singlet@meta.data$active.idents <- singlet@active.ident

HM <- data.frame(singlet[[colnames(ES)]], Idents(singlet)) ; colnames(HM)[ncol(HM)] <- condition
HM.Result <- getSignificance(HM, group = condition, fit = "Wilcoxon") 
HM.Result <- HM.Result[order(HM.Result$FDR),]
HM.Result <- HM.Result[HM.Result$FDR<0.05,]

HM.SC.Result.raw.RCHOP <- HM.Result
save(HM.SC.Result.raw.RCHOP, file = "../Reponse/HM.SC.Result.raw.RCHOP.Rdata")


HM.RCHOP.Result <- SC.Enrichissement(singlet, ES, condition)
DT.table(HM.RCHOP.Result)

Excipient

load(file = "../Reponse/Excipient_ES.RDs")

singlet.exp <- AddMetaData(singlet.exp, ES)
singlet.exp@meta.data$active.idents <- singlet.exp@active.ident

HM <- data.frame(singlet.exp[[colnames(ES)]], Idents(singlet.exp)) ; colnames(HM)[ncol(HM)] <- condition
HM.Result <- getSignificance(HM, group = condition, fit = "Wilcoxon") 
HM.Result <- HM.Result[order(HM.Result$FDR),]
HM.Result <- HM.Result[HM.Result$FDR<0.05,]

HM.SC.Result.raw.Excipient <- HM.Result
save(HM.SC.Result.raw.Excipient, file = "../Reponse/HM.SC.Result.raw.Excipient.Rdata")


HM.Excipient.Result <- SC.Enrichissement(singlet.exp, ES, condition)
DT.table(HM.Excipient.Result)

Pré-greffe

load(file = "../Reponse/Pre-greffe_ES.RDs")

singlet.pg <- AddMetaData(singlet.pg, ES)
singlet.pg@meta.data$active.idents <- singlet.pg@active.ident

HM <- data.frame(singlet.pg[[colnames(ES)]], Idents(singlet.pg)) ; colnames(HM)[ncol(HM)] <- condition
HM.Result <- getSignificance(HM, group = condition, fit = "Wilcoxon") 
HM.Result <- HM.Result[order(HM.Result$FDR),]
HM.Result <- HM.Result[HM.Result$FDR<0.05,]

HM.SC.Result.raw.PG <- HM.Result
save(HM.SC.Result.raw.PG, file = "../Reponse/HM.SC.Result.raw.PG.Rdata")

HM.PG.Result <- SC.Enrichissement(singlet.pg, ES, condition)
DT.table(HM.PG.Result)

Intersection

Diagram

Liste

Reduce(intersect, list(HM.PG.Result$Hallmark,HM.RCHOP.Result$Hallmark,HM.Excipient.Result$Hallmark))
##  [1] "HALLMARK_MITOTIC_SPINDLE"                  
##  [2] "HALLMARK_TNFA_SIGNALING_VIA_NFKB"          
##  [3] "HALLMARK_G2M_CHECKPOINT"                   
##  [4] "HALLMARK_HEME_METABOLISM"                  
##  [5] "HALLMARK_INTERFERON_GAMMA_RESPONSE"        
##  [6] "HALLMARK_MYOGENESIS"                       
##  [7] "HALLMARK_MTORC1_SIGNALING"                 
##  [8] "HALLMARK_UNFOLDED_PROTEIN_RESPONSE"        
##  [9] "HALLMARK_COAGULATION"                      
## [10] "HALLMARK_ESTROGEN_RESPONSE_EARLY"          
## [11] "HALLMARK_APICAL_JUNCTION"                  
## [12] "HALLMARK_IL2_STAT5_SIGNALING"              
## [13] "HALLMARK_KRAS_SIGNALING_DN"                
## [14] "HALLMARK_IL6_JAK_STAT3_SIGNALING"          
## [15] "HALLMARK_KRAS_SIGNALING_UP"                
## [16] "HALLMARK_HEDGEHOG_SIGNALING"               
## [17] "HALLMARK_HYPOXIA"                          
## [18] "HALLMARK_PROTEIN_SECRETION"                
## [19] "HALLMARK_APOPTOSIS"                        
## [20] "HALLMARK_COMPLEMENT"                       
## [21] "HALLMARK_ANDROGEN_RESPONSE"                
## [22] "HALLMARK_CHOLESTEROL_HOMEOSTASIS"          
## [23] "HALLMARK_ALLOGRAFT_REJECTION"              
## [24] "HALLMARK_PI3K_AKT_MTOR_SIGNALING"          
## [25] "HALLMARK_ANGIOGENESIS"                     
## [26] "HALLMARK_INFLAMMATORY_RESPONSE"            
## [27] "HALLMARK_PEROXISOME"                       
## [28] "HALLMARK_XENOBIOTIC_METABOLISM"            
## [29] "HALLMARK_E2F_TARGETS"                      
## [30] "HALLMARK_MYC_TARGETS_V2"                   
## [31] "HALLMARK_MYC_TARGETS_V1"                   
## [32] "HALLMARK_INTERFERON_ALPHA_RESPONSE"        
## [33] "HALLMARK_UV_RESPONSE_UP"                   
## [34] "HALLMARK_GLYCOLYSIS"                       
## [35] "HALLMARK_REACTIVE_OXYGEN_SPECIES_PATHWAY"  
## [36] "HALLMARK_DNA_REPAIR"                       
## [37] "HALLMARK_EPITHELIAL_MESENCHYMAL_TRANSITION"
## [38] "HALLMARK_OXIDATIVE_PHOSPHORYLATION"        
## [39] "HALLMARK_FATTY_ACID_METABOLISM"            
## [40] "HALLMARK_BILE_ACID_METABOLISM"             
## [41] "HALLMARK_ADIPOGENESIS"

3. Figures

RCHOP

PCA

Seurat::DimPlot(object = singlet, group.by = condition, label.size = 5, pt.size = 1,reduction = "pca", label = T) & Seurat::NoLegend() & 
  xlab(label = paste0("PCA 1 : ", round(Seurat::Stdev(singlet[["pca"]])[1],2), " %")) & ylab(label = paste0("PCA 2 : ", round(Seurat::Stdev(singlet[["pca"]])[2],2), " %"))

Volcano plot

SC.Volcano(response, FC)

Heatmap

Reponse
DoHeatmap(singlet, features = response$gene[1:50], group.by = condition, assay = "RNA")

Sample
DoHeatmap(singlet, features = response$gene[1:50], group.by = "Sample", assay = "RNA")

Excipient

PCA

Seurat::DimPlot(object = singlet.exp, group.by = condition, label.size = 5, pt.size = 1,reduction = "pca", label = T) & Seurat::NoLegend() & 
  xlab(label = paste0("PCA 1 : ", round(Seurat::Stdev(singlet.exp[["pca"]])[1],2), " %")) & ylab(label = paste0("PCA 2 : ", round(Seurat::Stdev(singlet.exp[["pca"]])[2],2), " %"))

Volcano plot

SC.Volcano(response.exp, FC)

Heatmap

Reponse
DoHeatmap(singlet.exp, features = response.exp$gene[1:50], group.by = condition, assay = "RNA")

Sample
DoHeatmap(singlet.exp, features = response.exp$gene[1:50], group.by = "Sample", assay = "RNA")

Pré-greffe

PCA

Seurat::DimPlot(object = singlet.pg, group.by = condition, label.size = 5, pt.size = 1,reduction = "pca", label = T) & Seurat::NoLegend() & 
  xlab(label = paste0("PCA 1 : ", round(Seurat::Stdev(singlet.pg[["pca"]])[1],2), " %")) & ylab(label = paste0("PCA 2 : ", round(Seurat::Stdev(singlet.pg[["pca"]])[2],2), " %"))

Volcano plot

SC.Volcano(response.pg, FC)

Heatmap

Reponse
DoHeatmap(singlet.pg, features = response.pg$gene[1:50], group.by = condition, assay = "RNA")

Sample
DoHeatmap(singlet.pg, features = response.pg$gene[1:50], group.by = "Sample", assay = "RNA")

4. Genes up-régulés

RCHOP

RP

RP.up.RCHOP <- FindMarkers(singlet, ident.1 = item1, ident.2 = item2 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod", only.pos = T)
RP.up.RCHOP <- SC.Result(RP.up.RCHOP)
DT.table(RP.up.RCHOP)

RC

RC.up.RCHOP <- FindMarkers(singlet, ident.1 = item2, ident.2 = item1 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod", only.pos = T)
RC.up.RCHOP <- SC.Result(RC.up.RCHOP)
DT.table(RC.up.RCHOP)
rm(singlet) ; gc(); gc(); gc(); gc(); 
##            used  (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.13e+07   601   1.72e+07   920 1.72e+07   920
## Vcells 1.38e+09 10508   3.04e+09 23186 2.97e+09 22636
##            used  (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.13e+07   601   1.72e+07   920 1.72e+07   920
## Vcells 1.38e+09 10508   3.04e+09 23186 2.97e+09 22636
##            used  (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.13e+07   601   1.72e+07   920 1.72e+07   920
## Vcells 1.38e+09 10508   3.04e+09 23186 2.97e+09 22636
##            used  (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.13e+07   601   1.72e+07   920 1.72e+07   920
## Vcells 1.38e+09 10508   3.04e+09 23186 2.97e+09 22636

Excipient

RP

RP.up.Excipient <- FindMarkers(singlet.exp, ident.1 = item1, ident.2 = item2 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod", only.pos = T)
RP.up.Excipient <- SC.Result(RP.up.Excipient)
DT.table(RP.up.Excipient)

RC

RC.up.Excipient <- FindMarkers(singlet.exp, ident.1 = item2, ident.2 = item1 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod", only.pos = T)
RC.up.Excipient <- SC.Result(RC.up.Excipient)
DT.table(RC.up.Excipient)
rm(singlet.exp) ; gc(); gc(); gc(); gc(); 
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636

Pré-greffe

RP

RP.up.PG <- FindMarkers(singlet.pg, ident.1 = item1, ident.2 = item2 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod", only.pos = T)
RP.up.PG <- SC.Result(RP.up.PG)
DT.table(RP.up.PG)

RC

RC.up.PG <- FindMarkers(singlet.pg, ident.1 = item2, ident.2 = item1 , logfc.threshold = log2(FC), assay = "RNA", slot = "data", test.use = "bimod", only.pos = T)
RC.up.PG <- SC.Result(RC.up.PG)
DT.table(RC.up.PG)
rm(singlet.pg) ; gc(); gc(); gc(); gc(); 
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636
##            used (Mb) gc trigger  (Mb) max used  (Mb)
## Ncells 1.12e+07  599   1.72e+07   920 1.72e+07   920
## Vcells 8.67e+08 6617   2.43e+09 18549 2.97e+09 22636

5. Tableaux de synthèse

synthese <- data.frame(row.names = c('Excipient', 'RCHOP', 'Pré-Greffe'),
  DEG = c(length(response.exp$gene), length(response$gene),length(response.pg$gene)),
  UP = c(length(RC.up.Excipient$gene), length(RC.up.RCHOP$gene), length(RC.up.PG$gene)),
  DOWN = c(length(RP.up.Excipient$gene), length(RP.up.RCHOP$gene), length(RP.up.PG$gene))
)

DT::datatable(t(synthese), class = 'cell-border stripe')